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Abstract: We report on the response of a high light-output Nal(Tl) crystal to nuclear 
recoils induced by neutrons from an Am-Be source and compare the results with the response 
to electron recoils produced by Compton-scattered 662 keV 7 -rays from a ^^^Cs source. The 
measured pulse-shape discrimination (PSD) power of the Nal(Tl) crystal is found to be 
significantly improved because of the high light output of the Nal(Tl) detector. We quantify 
the PSD power with a quality factor and estimate the sensitivity to the interaction rate for 
weakly interacting massive particles (WIMPs) with nucleons, and the result is compared 
with the annual modulation amplitude observed by the DAMA/LIBRA experiment. The 
sensitivity to spin-independent WIMP-nucleon interactions based on 100 kg-year of data 
from Nal detectors is estimated with simulated experiments, using the standard halo model. 
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1 Introduction 

The existence of non-baryonic cold dark matter has been widely supported by many astro¬ 
nomical observations [1-6]. The theoretically favored dark matter candidates are weakly 
interacting massive particles (WIMPs), which are well motivated by supersymmetric mod¬ 
els [7, 8]. In the constrained minimal supersymmetric standard model, the lightest super- 
symmetric particle is a WIMP candidate (y) with an expected mass of > 100 GeV/c^ [9]. 
Recently, a number of experimental observations [10-13] could be interpreted as direct de¬ 
tection signatures of WIMP-nucleon interactions. Moreover recent astronomical gamma-ray 
observations could be interpreted as indirect signatures of WIMP annihilations near the cen¬ 
ter of the Galaxy [14]. Therefore unambiguous tests of the direct WIMP detection signals 
are needed. 

Among the direct searches for WIMPs, the DAM A/LIBRA experiment is unique in that 
it has been consistently reporting positive signals for an annual modulation of the nuclear 
recoil detection rate since 1998 [15], most recently with a 9.3cr statistical significance [10, 16, 
17]. This annual modulation is consistent with WIMP scattering off nuclei inside Nal(Tl) 
detectors. If the annual modulation signal observed by DAMA/LIBRA is interpreted in 
the context of the standard halo model [18-20], the favored regions of WIMP-nucleon 
cross-section and WIMP-mass parameter space have been excluded by other direct search 
experiments [21-26]. This has motivated alternative explanations for the DAMA/LIBRA 
signal [27-30]. However, when non-trivial systematic differences in detector responses [31- 
34] and reasonable modifications to the commonly used astronomical model for the WIMP 
distribution [35, 36] are considered, it remains possible to reconcile both the DAMA signal 
and the null observations [37, 38]. Therefore, it is important to devise experiments that 
search for the DAMA/LIBRA signal using the same technique but with higher sensitivity. 
To date, no independent, direct experimental verification of the DAMA/LIBRA signal has 
been performed. 

To reproduce the DAMA/LIBRA observation, several experimental groups have started 
developing ultra-low-background Nal(Tl) crystals [39-41] with target background levels at 
or below those measured in the DAMA/LIBRA experiment. If this background reduction is 
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successful, sizable experiments using more than 100 kg of crystals will be performed within 
the next few years. Accumulating a similar amount of data as that acquired with the 
DAM A/LIBRA experiment can test the DAMA/LIBRA annual modulation observation. 
However, this will require a long-term stable data taking over a period of at least three 
years. 

On the other hand, we can extract a WIMP signature based on differences between 
the scintillation characteristics of nuclear and electron recoils using a so-called pulse-shape 
discrimination (PSD) analysis. This technique has been applied for WIMP searches with 
CsI(Tl) crystals [25, 42, 43] and Nal(Tl) crystals [44, 45]. In the past, the discrimination 
power of Nal(Tl) crystals has been much poorer than that of CsI(Tl) crystals [46-48]. 
However recently developed Nal(Tl) crystals grown by Alpha Spectra Inc. provide higher 
light yields [39, 41], and this makes it interesting to study their PSD performance and 
evaluate their sensitivity. 

In this article, we present a new measurement of the PSD capability of a high light- 
output Nal(Tl) crystal. We measure the response to nuclear recoils induced by neutrons 
from a 300-mCi Am-Be neutron source and compare the results with the response to electron 
recoils induced by a ^^^Cs 7 source. The Seoul National University neutron calibration 
facility [48] is used for this measurement. Using the results from these measurements, 
we estimate the sensitivity of a Nal(Tl) dark matter search experiment that uses PSD. A 
model-independent comparison between the expected upper limit on the nuclear recoil rate 
from the new Nal(Tl) experiment and the DAM A annual modulation amplitude promises 
a better understanding of the underlying scenarios of current dark matter search results. 

2 Measurement of the pulse-shape discrimination 

For the neutron-nuclei elastic scattering measurements, the neutron source was surrounded 
by liquid scintillator (LSC), lead, and polyethylene shields with a collimated aperture in the 
crystal direction as shown in Fig. 1. A small 2x2xl.5 cm^ Nal(Tl) crystal made from the 
same. Alpha Spectra-grown, ingot as a large crystal (Nal 002) used for background mea¬ 
surements in the Yangyang underground laboratory [41] was used as a target. Two 3-inch 
photomultiplier tubes (PMTs) that were selected for high quantum efficiency (RI2669SEL, 
Hamamatsu Photonics) were attached to each end of the Nal(Tl) target crystal, which was 
located at the exit of the collimated aperture of the neutron source. To identify neutrons 
that elastically scatter in the crystal, two neutron-tagging detectors, consisting of BC50IA 
liquid scintillator contained in cylindrical 0.5 i stainless-steel containers, were located at 
90° from the line between the neutron source and the Nal(Tl) crystal. The amplified signals 
from the crystal and the neutron detectors were recorded by 400-MHz flash analog-to-digital 
converters for 25-/is time intervals. The trigger condition for the Nal(Tl) crystal was one 
or more photoelectrons in each PMT within a 200 ns window, which is the same condition 
that was used for the underground data described in Ref. [41]. We additionally required 
a time coincidence between the Nal(Tl) crystal and one of the neutron detectors to occur 
within a 2 /rs time window. For the electron recoil measurements, the crystal was irradiated 
with 662 keV 7 -rays from a ^^^Cs source. The data were taken separately with the same 
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Figure 1. A setup for the neutron-nuclei elastic scattering measurements using the Am-Be neutron 
source. 


setup, except for a lead and polyethylene beam stopper that was placed in the collimated 
aperture of the neutron source. The trigger condition for the Nal(Tl) crystal was the same; 
the neutron detectors were not used. 

The energy calibration of the Nal(Tl) crystal was performed with 59.54 keV 7 -rays 
from an ^^^Am source. A clustering algorithm described in Ref.s [42, 48] is applied to 
separate single photoelectrons as individual clusters and reduce the low energy baseline 
noise. The measured electron-equivalent energy was determined from the charge sum of the 
identified clusters, and we obtain a light yield for this crystal of approximately 16 photo- 
electrons/keV, which is consistent with the measured yield for large crystals grown by the 
same company [39, 41]. This light yield is approximately twice larger than that of crystals 
that were used in previous studies [46, 49, 50]. 

In the offline analysis, we require a tighter coincidence condition of less than approx¬ 
imately 100 ns between the signals from the Nal(Tl) crystal and the neutron detectors, 
as discussed in Ref. [48]. Neutron-induced signals in the neutron detectors are identified 
by taking advantage of the nuclear and electron recoil PSD capabilities of BC501A liquid 
scintillator [51, 52]. We recorded neutron data with the Am-Be neutron source over a 
three-month period. 

Because of the different time distributions of the photoelectrons produced by nuclear 
and electron recoils in scintillation crystals, it is possible to use PSD to distinguish the 
nuclear recoil signals from the electron-recoil-dominant backgrounds. To characterize the 
time distributions, we use the natural logarithm of the mean decay time (In(MT)), defined 
as 

In(MT) = In - ^o) , (2.1) 

where Ai is the charge of the ith cluster, ti is the time of the ith cluster, and to is the time 
of the first cluster. Here, only clusters within 1.5/rs from to are considered. Figure 2 shows 
the In(MT) distributions for the nuclear and electron recoil events in the Nal(Tl) crystal 
for 1 keV bins of visible energy for the range 1 < E < 5 keV. In this figure, differences 
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(a) 1-2 keV 



(c) 3-4 keV 



(b) 2-3 keV 



(d) 4-5 keV 


Figure 2. The In(MT) distribution of neutron-induced events (open circle) from the Am-Be 
source and electron-induced events (filled circles) from the source. The results of fits using 

asymmetric Gaussian functions are superimposed. 


in the In(MT) distributions between nuclear recoils from electron recoils are evident, even 
for energies as low as 1 keV. This result is significantly improved compared with those of 
previous studies that reported PSD capabilities only for energies of 4 keV and higher [44- 
46]. This gives us confidence to design an experiment for WIMP searches with Nal(Tl) 
crystals based on the PSD performance. 

We performed fits to the In(MT) distributions using asymmetric Gaussian functions 
that are shown in the plots, where good agreement with the data is observed. Similar fits 
were performed for each 1 keV energy bin from 1 to 14 keV for nuclear and electron recoil 
events, and Fig. 3 shows the most probable value and the root mean square (RMS) of 
In(MT) for electron and nuclear recoil data. As shown in the plots, the separation power is 
reduced at low energy, but some PSD capability persists at energies as low as 1 keV, which 
is even better than CsI(Tl) crystals [47, 48]. 

We characterize the PSD power using a quality factor [53] that is commonly used to 
characterize WIMP-search detectors [46-48] and defined as 
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Figure 3. (a) The most probable value of In(MT) from the asymmetric Gaussian fits for each en¬ 

ergy bin for neutron-induced (open circles) and electron-induced (filled circles) events, respectively, 
(b) The root mean square (RMS) of In(MT) as a function of energy for neutron-induced (open 
circle) and electron-induced (hlled circle) events. 


where a and /3 are fractions of the nuclear and electron recoil events that satisfy the selection 
criteria, respectively, where we have varied these criteria to provide the best (he., the 
minimum) value of the quality factor. For a detector with perfect discrimination, a = 1 
and /3 = 0, thus a smaller quality factor means a better PSD power. Figure 4 shows 
the measured quality factor of the Nal(Tl) crystal studied here compared with results 
from previous Nal(Tl) [46] and CsI(Tl) crystal measurements [47, 48]. As shown in this 
figure, the quality factor of this new, high light-yield Nal(Tl) crystal is approximately one 
order of magnitude lower than that from previous studies [46]. To demonstrate that this 
is mainly due to the approximately two times higher light yield of the new crystal, we 
adjust our results by assuming a light yield similar to typical Nal(Tl) crystals and find 
results that agree with previous measurements. We also find that the quality factor of this 
crystal is slightly better than that of CsI(Tl) crystals, especially at energies below 6 keV. 
This promises relatively good sensitivity for low-mass WIMP detection with these Nal(Tl) 
crystals. 

3 Expected sensitivity of a Nal(Tl) experiment using PSD 

The Korea Invisible Mass Search (KIMS) Collaboration recently joined a consortium that is 
developing low-background, high light-yield Nal(Tl) crystals for attempts to reproduce the 
DAM A/LIBRA observation. We have achieved a background level of approximately 3 dru at 
6 keV with a below 2 keV energy threshold [41]. Because we already identified the significant 
backgrounds of the Nal(Tl) crystals as being due to ^^'^Pb and '^^K contaminations, the 
background levels of future crystals that are produced as part of this program will be 
significantly reduced. Even without internal background reduction, it should be possible 
to achieve a ~2 dru background level by immersing the Nal(Tl) crystal array inside a 
liquid scintillator box that serves as an active veto. After growing Nal crystals with a 
total mass of 100 kg, which may be completed by the end of 2015, the KIMS collaboration 
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Figure 4. The measured quality factor of the new Nal(Tl) crystal compared to previous Nal(Tl) 
measurement by Gerbier et al. [46] and CsI(Tl) measurements using monoenergetic neutron 
beam [47] as well as the Am-Be source [48] . 


plans to start the KIMS-Nal experiment that, after three years of stable operation, could 
unambiguously confirm or refute the DAMA/LIBRA annual modulation signals without 
any model dependence. 

In the meantime, we plan to extract information about the WIMP interaction using 
Nal(Tl) crystals as early as possible. One of our goals is to extract the WIMP interaction 
signal using data from the initial period of one year or less. Because of the very good PSD 
power of the new Nal(Tl) crystals, the extraction of a meaningful cross section limit for 
WIMP-nucleon interactions in Nal should be possible with a small data set. To quantify the 
performance of this measurement, we conservatively estimate the sensitivity of the KIMS- 
Nal experiment by assuming a 100 kg Nal(Tl) crystal array with a 2 dru background, a 
1-year exposure, and a 2 keV energy threshold. 

We performed simulated experiments to extract the upper limits of the nuclear recoil 
events. For events in each 1 keV bin, we ran 2,000 simulated experiments based on the 
expected electron recoil background, assuming no nuclear recoil signal events. We construct 
Bayesian likelihoods that are formed as products of the measured log(MT) probability 
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Figure 5. (a) The median expected 90%-confidence-level upper limits (points) with Icr standard 

deviation probabilities (error bars) on the rate of nuclear recoil events from the KIMS-Nal 1- 
year experiment (filled circle) are compared with the annual modulation amplitudes observed by 
the DAM A/LIBRA experiment (open circle), (b) Expected 90%-conhdence-level upper limit on 
the annual modulation amplitude as a function of the annual modulation ratio of WIMP-nucleon 
interaction from the KIMS-Nal experiment, which is compared with the DAMA/LIBRA observation 
in the 2-6 keV energy range. 


distributions for nuclear recoil signals and electron recoil backgrounds. For each simulated 
experiment, 90% confidence level (CL) limits are determined such that 90% of the posterior 
densities of the nuclear recoil event rates fall below the limit. In Fig. 5 (a), we present 
the median expected 90% CL upper limits (points) together with Icr standard deviation 
probability (error bars) of the nuclear recoil events for each 1 keV bin. The DAMA/LIBRA 
annual modulation amplitudes for each 0.5 keV bin from Ref. [19] are overlaid in the plot 
for comparison. At low energies, the KIMS-Nal expected upper limits of the nuclear recoil 
rates corresponding to the WIMP interaction with Nal(Tl) nuclei are consistent with the 
annual modulation amplitudes of the DAMA/LIBRA experiment. 

To make a direct comparison between the expected nuclear recoil rate from the KIMS- 
Nal experiment and the annual modulation amplitudes observed by the DAMA/LIBRA 
experiment, we need to know the expected ratio of the annual modulation amplitude to the 
total WIMP-nucleon interaction rate, which depends on the models of the halo structure 
as well as the underlying theory of the WIMP candidate. In other words, a WIMP-nucleon 
interaction limit determined using the PSD analysis, can constrain the ratio of the an¬ 
nual modulation amplitude to the total rate of the WIMP-nucleon interaction, which is 
model-independent. The annual modulation amplitude observed by the DAMA/LIBRA 
experiment for the 2-6 keV events with the highest significance is 0.0112±0.0012 dru [10]. 
The expected 90% CL upper limit of the nuclear recoil rate from the KIMS-Nal experiment 
is 0.0104±0.0050 dru. In Fig. 5 (b), we convert the expected upper limits of the nuclear 
recoil rate from the KIMS-Nal experiment into an annual modulation amplitude as a func¬ 
tion of the annual modulation ratio from the total WIMP-nucleon interaction rate. We find 
that the KIMS-Nal experiment can investigate the allowed parameter space of the DAMA 
2cj signal region for any model with an annual modulation ratio less than 85%, assuming 
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Figure 6. The median expected (black dotted line) 90%-confidence-level upper limit on the 
WIMP-nucleon spin-independent cross-section assuming the background-only hypothesis, shown 
together with WIMP-induced DAM A/LIBRA 3cr-allowed region (solid contour). The dark (green) 
and light (yellow) shaded bands indicate the Icr and 2a standard deviation probability regions 
in which the limits are fluctuated. The red dotted line indicates the expected limit assuming an 
optimistic background level of 1 dru and an energy threshold of 1 keV. 


the median expectation shown in Fig. 5 (b). 

Assuming the standard halo model [18], we estimate the sensitivity of the KIMS-Nal 
experiment as the predicted cross-section limits for the WIMP-nucleon spin-independent 
interactions in case of no signal for an ensemble of Monte-Carlo experiments. For each 
experiment, we determine simulated log(MT) distribution for an electron-recoil-background 
only hypothesis. The rate of background, assuming a 1-year exposure of a 100 kg array 
of Nal(Tl) detectors with 2 dru background rate, has Poissonian fluctuations. We then fit 
the log(MT) spectrum with signal plus background hypothesis with flat priors for both the 
signal and the background rate. We construct a Bayesian likelihood for each 1 keV bin 
from fits to the log(MT) distribution. The overall likelihood is obtained by multiplying 
the likelihood for each bin between 2 and 10 keV considering the energy-dependent rate of 
WIMP interaction for WIMP masses between 5 GeV/and 10000 GeV/. The 90% CL 
limit for each simulated experiment is determined such that 90% of the posterior density 
of the WIMP-nucleon cross-section falls below the limit. The median expected 90% CL 
limits and la and 2a standard-deviation probability regions are calculated from the result 
of 2,000 simulated experiments. To directly compare these results with the DAMA/LIBRA 
allowed signal region, we use quenching factors from Ref. [19] of 0.3 and 0.09 for Na and 
I, respectively. The quenching factor is the fraction of the recoil energy deposited by a 










WIMP with respect to the energy deposited by an electron (or 7 -ray). Figure 6 shows the 
expected median limit, and Icr and ‘la standard deviation probability regions, of the KIMS- 
Nal experiment using PSD analysis. We also include results from an optimistic KIMS-Nal 
expected limit of a 1 dru background rate with a 1 keV energy threshold. These results 
are further compared with the DAMA/LIBRA 3a allowed signal region from Ref. [19]. 
As shown in this figure, the median expected limit from the KIMS-Nal 100 kg-year data 
can investigate the DAMA/LIBRA 3a regions even for the conservative threshold and 
background level assumptions. This comparison, based on measurements using the same 
target-detector material, may provide a better understanding of the nature of dark matter 
and of the DAMA/LIBRA signature. If the KIMS-Nal experiment obtains similar annual 
modulation signatures using Nal(Tl) crystals in the future, this PSD measurement will 
provide valuable information about the annual modulation events that will help determine 
whether or not they originate from WIMP-induced nuclear recoil events. 

4 Conclusion 

We report measurements of the PSD characteristics for the nuclear and electron recoil events 
of the Nal(Tl) crystal. Because of the high light output of recently developed Nal(Tl) 
crystals, we find a PSD quality factor that is approximately one order of magnitude better 
than those of previously studied crystals. The influence of this good PSD capability on 
the expected sensitivity of the KIMS-Nal experiment is evaluated. Model independent and 
dependent comparisons of the DAMA/LIBRA annual modulation signature give confidence 
that the KIMS-Nal experiment can investigate interesting regions of parameter space with 
one year of accumulated data. 
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